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TECHNICAL NOTE 3507 

PRACTICAL CONSIDERATIONS IN SPECIFIC APPLICATIONS 
OF GAS-FLOW INTERFEROMETRY 
By Walton L. Howes and Donald R. Buchele 


SUMMARY 

Optical refraction, end and comer effects, and spurious inter- 
ferences may significantly affect interferometric evaluation of gas- 
density fields. Evaluation equations wMcH account for refraction are 
derived in a previous report. In the present report these equations 
are extended to simultaneously include corner effects. The spurious 
interferences can be used to aline a model with the light beam. A 
simple method for evaliiating the density adjacent to a surface is 
described. 

The random error in measuring fringe shifts is a fimction of fringe 
spacing. In a representative experiment this random error was found to 
be less than the fringe- spacing variation produced by free-stream t\ir- 
bilLence and optical imperfections. The latter variation was consider- 
ably less than the fringe shifts caused by steady- state density varia- 
tions associated with boundary layers. 

Recomputed la m inar-bovindary- layer density and velocity profiles 
associated with supersonic flow along an Insulated flat plate are in 
fairly good agreement with the theory of Chapman and Rubesin. Computed 
skin-friction coefficients are in fair agreement with theory. The re- 
sidual disagreement between theory and experiment is attributed to 
thickening of the boundary layer in the vicinity of the plate midspan 
plane. Temperature variations within the windows may also be significant. 

Interferometric determinations of the density at the wall of a flat 
plate in subsonic flow are in very good agreement with results given by 
other methods . 

Interferometric determinations of free- convection temperature pro- 
files beneath a heated horizontal cylinder are in very good agreement 
with the theory of Hermann when the kinematic viscosity of air is based 
on the cylinder wall temperature, and refraction and corner effects are 
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accounted for. Agreement is poor, especially for large temperature dif- 
ferences, when the kinematic viscosity is based upon ambient temperature. 


According to the eAjjcrxiucuual results, refraction is generally 
significant in boundary-layer studies. Corner effects are likely to 
be important when the model is not bounded by windows, but unimportant 
when the model is bounded. In determining the free-stream density in 
a wind tunnel, end effects are usually significant. 

In three typical experiments the applicability, apparent-ray- 
trace crossing, and source- size criteria specified in a previous report 
were found to be satisfied. Series remainders were generally negligible. 


Density evaluation equations for nse with optical interferograms 
of one-dimensional density fields are derived in reference 1. These 
equations accoiont for optical refraction. However, certain additional 
sources of error are of interest. The most significant additional 
sources of error are likely to be 

(1) End and corner effects 

(2) Measurement of interference-fringe shifts 

(3) Model allnement 

(4) Spmrious fringe shifts in the vicinity of surfaces 

(5) Turbulence 

All are considered in the present report. 

Refraction corrections are described in references 1 to 3. Correc- 
tions for end or corner effects are reported in references 4 to 9. 
However, optical refraction and corner effects have not been analyzed 
simultaneously, although the two effects are interdependent. In the 
present report the evaluation equations reported in reference 1 are ex- 
tended to Include simultaneous corrections for corner effects. 

In certain instances, the usefulness of interferometry may depend 
upon the acciiracy of interference-fringe-shift measurements. Therefore, 
the resvilts of an investigation of the random errors in measuring fringe 
shifts are included in the present report. 
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In addition^ spurious fringe shifts may introduce significant 
errors in evaluating the density immediately adjacent to a surface. 
Sources and uses of the spurious fringe shifts are considered. A tech- 
nique not involving spurious interferences is described for evaluating 
the density adjacent to a surface. The effect of turbulence is noted. 

The significance of the preceding items, as well as other consid- 
erations indicated in reference 1, is determined for three representa- 
^ tive experimental studies, namely, 

(1) The laminar boundary layer associated with supersonic air flow 
along an insulated flat plate 

(2) The boimdary- layer density distribution associated with sub- 
sonic air flow along a thick flat plate with a blunt leading 
edge 

(3) The temperature distribution associated with free- convective 
heat transfer from a horizontal cylinder 


MEASUREMENT OF FRINGE SHIFTS 

Suppose that in a given experiment the density field recorded as 
an optical interf erogram consists of an ambient region having constant 
density p ^ and a second region in which the density p is a fimction 
of a single Cartesian coordinate y (and time) which is essentially 
perpendicular to the path of the interferometer light beam. The inter- 
ference fringes are, then, usxaally initially oriented 

(1) Parallel to the gradient of p 

(2) Perpendicular to the gradient of p 
or 

(3) For the infinite-fringe condition 

Orientations 1 and 2 are iJlLustrated in figures l(a) and (b) , respec- 
tively. Orientation 3 corresponds to constant phase of interference. 

The fringe pattern associated with the field p(y) differs from 
that which would exist if p(y) were replaced by p^ (compare corre- 
sponding "flow" and "no flow" patterns in fig. l). At any point on 
the interferogram the change of phase and order of interference cor- 
responding to the change p^ to p(y) is called the fringe shift at 
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that point. Thus^ fringe shifts^ henceforth denoted by are ex- 
pressed in units of fringes^ where a unit fringe shift corresponds to 
a change of 1 order of interference. (All symbols are defined in 
appendix A.) 

Expressions for computing the fringe shifts from measurements of 
lengths on an interferogram are as follows: 

(l) Orientation parallel to gradient of p - At a given value of 
y^ (where primed quantities^ henceforth^ refer to lengths on the inter- 
ferogram and^ thus^ ta^e into account magnification)^ 

/o < |x' - XjJil < 

. n . in, 0,1,2,.. 

d' \ 

~ \ n < I m 

where x' - x^ (in the image plane) is the displacement in the x' di- 
rection of a fringe with respect to the position Xq of a given fringe 
for ambient conditions, and d^ is the mean fringe pitch associated 
with the field p . A simple accurate method for determining any value 

oo 

of N is suggested by the preceding formula and will be described sub- 
sequently. The random error AN in N is given by 





(N / 0) 


where Ad^ is the random error in a single measurement of d^ . 

oo OD 

(2) Orientation perpendicular to gradient of p - At a given 
value of y ' , 


iNl = n - i ^d; (v = 1, 2, 3, . . n) 

^«o V=1 

where the d(, are associated with the field p(y) and are s\ammed 
starting from the value of y at which p(y) = p^. The quantities 
d' and d’ may represent either fringe pitch or semipitch. Only 
V * 1 1 

values d^, d^ = 1> • • • are included, because only interfer- 

ence maximums and minimums can be determined reliably by visual observa- 
tion. The random error in N is given by 
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(N 4 0) 


(3) Infinite-fringe adjustment - At a given value of y' , 


|W| = n 


where 



The random error in N is given hy 



(K i 0) 


where the error of estimating the phase of the infinite fringe associ- 
ated with is estimated to be about l/S fringe, and Ad^ corresponds 

to the value of Ad' when d' = d' . 


The question arises as to which initial fringe adjustment is most 
desirable. Initial orientation parallel to the gradient of p appears 
to be preferable for the following reasons: 

(1) Accurately measurable values of N are not limited to spe- 
cific values d' = 1, 1^, ..... r 

(2) It follows from reason 1 that p(y) can be determined rela- 

tively accin-ately even if the maximum absolute value of N amoimts to 
less than 1 fringe. ^ ^ 

(3) The possibility of confusing an interference fringe adjacent 
to a sxirface with the surface itself is reduced. 

(4) The possibility of confusing interference fringes adjacent to 
a surface with diffraction fringes parallel to the s\rrface is reduced. 


Advantage 1 applies specifically to visual measurements on an in- 
terferogram. Photoelectric measurements on an interferogram also allow 
determination of intermediate values of N for the other orientations. 
However, accurate photoelectric measiirements require extremely precise 
photographic technique because photographic exposure is a nonlinear 
function of the phase of interference, and photographic density is 
generally a nonlinear function of photographic exposure. 
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For initial orientation parallel to the gradient of the pro- 
files of adjacent fringes are Identical. Thus^ for each selected 
value of y' the shift x' - of a given fringe from its ambient 

position is measured. The fringe shift is computed from the formula 
previously given^ where n = 0. However^ when x' - Xq attains a 
value such that |x' - x^| ^ d^^ the x' traverse of the comparator is 

returned (holding y' fixed) to the adjacent fringe, which has shifted 
to the vicinity of Xq. Measurements of x' - Xq for selected values 

of y' are made on this new fringe and N is computed with n = 1. 
When again |x' - x^j ^ d^ for this fringe, the next adjacent fringe 

is selected for measvirement and N is computed with n = 2. The 
process can be repeated for any number of fringe shifts. The measiired 
quantities are illustrated in figure l(a) . The process is extremely 
simple, and accurate measurements of arbitrary values of N are ob- 
tainable. All measurements are performed between the extensions of 
two adjacent ambient fringes. Thus, a fairly large x component of 
grad p can be tolerated. Intolerable values of the x component 
of grad p cause discontinuities of the measured fringe-shift profile 
at integer values of N. In such cases it would be necessary to re- 
vert to the method of measurement described in reference 10. 


RANDOM ERRORS IN DETERMINING FRINGE SHUTS AND FRINGE- SHUT DERIVATIVES 

The importance of random errors in measuring fringe shifts is dis- 
cussed in reference 11. The probable error in measuring fringe shifts 
was assumed to be 1 fringe (/IN = l-l) . The minimum required total fringe 
shift for accurate determination of a variable density field was esti- 
mated to be at least 10 fringes (N = i-lO) . However, more recent ad- 
vances cause the ass\mied random errors and conclusions reached in refer- 
ence 11 to appear extremely conservative . For example, the relative 
error in experimentally measuring N is reported to be i0.07 in refer- 
ence 8 and, from the data given, is about i0.04 to i0.09 in reference 
12. Consequently, the minim\jm reqviired total fringe shift is consider- 
ably less than 10 fringes, and the range of Mach numbers (in the case 
of flow phenomena) for which the interferometric method is useful is 
probably considerably greater than the range 0.5 to 3.0 predicted in ref- 
erence 11. 

The error in N is often regarded as the largest random error in 
gas-flow Interferometry. Moreover, according to the methods described 
in references 1 and 2, measiirements of slopes dN/dy^ of the measured 

fringe-shift profile may also be desirable. It becomes important to 
determine what the random errors in measxiring N and dN/dy^ actually 
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Error in N 

The random error in N associated with steady-state measurements 
of boundary layers and similar phenomena -is composed of two types of 
error ^ namely, 

(1) The random error involved in measuring N 

(2) The random error ca\ised by fluctuations of the phenomenon about 
the steady state 

The error in measuring N by visioal methods depends upon several 
factors, among which are 

(1) The value of d^ 

(2) The distribution p(y) 

(3) The initial fringe orientation with respect to the gradient 
of p 

(4) The photographic properties of the interference-recording 
medium, for example, contrast, density, resolution, graininess 

(5) The magnifying power of the instrument used for viewing the 
interferograms 

The effect of dj|^ is considered in the following paragraphs. The ef- 
fect of initial fringe orientation is presented in the preceding sec- 
tion. Items 2, 4, and 5 are discussed in appendix B. 

Consider two adjacent parallel interference fringes having a center- 
to-center separation (pitch) dj^. Then, the relative error in a single 

measurement of due to an error Ad^ in measuring d' is given by 


Ad' 

oo 



where^ as before^ Ad* represents the random error In a single measure- 

OO 

ment of d*^ d* represents the mean value of several measurements of 

“"i i“-i - 

The standard deviation in measuring N was determined from a 

OO 

series of 19 interferograms possessing various values of dj^ by means 
of a commercial optical comparator having the following characteristics: 
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(1) Readability; 0.001 millimeter (0.00004 in.) 

(2) Magnifying power: X22.6 

(3) Illumination; Diffused 

All measurements were made at one location on each interferogram in 
order to avoid the random error associated with spatial variation of 
d' . Additional details are given in appendix B. 

00 

The results of determining AdJ_^ and AN^ as functions of dj^ are 

shown in figures 2(a) and (bj[^ respectively. The limiting profile 
(ANjjiiiQ as a function of corresponding to the readability of the 

comparator^ is also shown in figure 2(b). The measurements of Ad^ 
are represented by the straight lines 


Ad^ « 0.003 mm (d* < 1.4 mm) 

Ad' “ O.OOSdLmm (d' >1.4 mm) 

00 OP 


so that 

AN^ » 0.003/d;, 

AN -0.005 

OP 

where Ad^ and d^ are in millimeters. The pair of profiles associ- 
ated with each of the above sets of formulas do not coincide at 
dj^ = 1.4 millimeters. This lack of coincidence is of minor importance^ 

however^ because only orders of magnitude are of interest. The minimum 
value * 0.005 is 1/200 of the value assumed in reference 11 and 

about l/lO of the values indicated in references 8 and 12. However, the 
values indicated in references 8 and 12 probably include all factors con- 
tributing to AN, except possibly time. 


(d^ <1.4 mm) 

OP 

(d' >1.4 mm) 

OP 


t 


Because of the smallness of Ad', it is quite likely that space 

OP 

and time variations in d;^ caused by turbulence and optical imperfec- 
tions will be greater than AdJ^ even over small regions of the field. 
Space variations of d^^ can be determined from a single interferogram. 

Time variations can be determined for corresponding points on a series 
of interf erograms . Space variations of d;^ were determined from single 

measurements at several locations on a single interferogram. The space 
variations (standard deviations) associated with free-stream region in 
the subsonic-flow experiment to be discussed were 
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Ad^ * 0.02 mm 
00 

AN *= 0.01 

oo 

for » 1.31 millimeters and within a fixed interval of x^. The 
latter variations are two or three times greater than the corresponding 
deviations associated with the measurement process. However^ the spa- 
tial variation AI^ « 0.01 is only about 1 percent of the maximum 
fringe shift involved in the experiments. Time variations in N are 
discussed in conjunction with the subsonic-flow experiment later in the 
present report. 

All results were obtained for j N^| = 1. If the interval of meas- 
urement includes n constant pitch fringes ^ then the relative error in 
measuring nN«o is inversely proportional to n^ because the error is 
influenced only by the observer's ability to determine the locus of the 
interference maximimis (or minimums) bounding the interval. 


Error in dw/dy^ 

The slope dN/dy^, of a measured fringe-shift profile Np, as re- 
quired by the evaluation procedure presented in reference 1, can be de- 
tennined by several different methods. The error in determining the 
slope depends upon the error in the profile Nj) as well as the error 

in measuring the slope. However, by assuming a known profile Np, tests 

were made to determine the error in measuring the slope by means of an 
optical differentiator. The optical differentiator was selected be- 
cause of its simplicity and the rapidity with which measurements can be 
performed. 

A differentiator based on the principle of that described in ref- 
erence 13 was constructed with provision for attachment to a drawing 
machine, thus permitting direct reading of slope angles co in incre- 
ments of 5 minutes of arc. Measurements of (a were made from a 7- by 
10-inch plot of the profile Njj shown in figure 3(a) . This profile 

was computed in reference 1 from a hypothetical exponential-density pro- 
iiis* The value of dW/dyj^ associated with any given abscissa value 
yj) is given by 


dH 


= (scale factor) tan to 


Comparison of the hypothetical and measured profiles dN^/dy^ is shown 
ill i*i§nre 3(b) . The standard deviation of a single measurement of cu 
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was found to vary from zero to about 10 minutes of arc^ the larger error 
being associated with the greatest c\rrvat\ire of the profile For 

the range of values of dU/dy^ considered, the relative error in de- 
termining dN/dyj) was less than 10 percent, except when 1 I approached 
jt/2, in which case small errors in measuring co correspond to very large 
errors in tan oi. 


Error in Density Profile 


With some knowledge available regarding random errors in measuring 
N and dN/dyj), the random errors in the density ratio 0^ and distor- 
tion D, which determine the density profile p(y), can be computed. 


Evaluation equations for p* and D are given in reference 1 and 
are repeated in appendix B. Expressions for the random errors in p* 
and D also are given in appendix B, where it is shown that the error 
in measuring N is of prime importance in determining the error in p* 
when I Nl is small but becomes much less significant as InI increases. 
Because the absolute magnitude of the refraction term, namely. 


(2 - 3K)b^L^, is much less than that of the primary term N i 

the eciviation for p*, it follows that relatively large errors in dw/dy^^ 


have only a small effect upon p*. For example, if the refraction term 
amoxmts to 10 percent of the primary term, then an error of 10 percent 
in determining dU/dy^, amounts to an error of less than 2 percent in 

computing p*. However, a 10-percent error in determining dH/dy^^ 


yields a 10-percent error in computing D. Thus, it is desirable to 
keep D as smal 1 as possible. This can be accomplished by properly 
choosing the object plane, that is, the value of K» Generally, the 
previously described method for determining dN/dy^ is adeq.uate. 


The preceding paragraphs deal with the errors in evaluating the 
shape of a profile p(y) without regard to its location with respect to 
material objects. Errors involved in locating the profile are somewhat 
different from those discussed previously and may best be considered in 
conjunction with specific experiments. 


SUEFACE EFFECTS AND MODEL ALINEMENT 

Quite often the region immediately adjacent to a model surface is 
of prime Interest in gas-flow experiments. Unfortunately, this region 
is likely to be the region in which experimental measurements of any 
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kind are most difficiolt. In the case of interferometry, the principal 
deleterious effects are misalinement of the model with respect to the 
light beam and spurious fringe shifts caused by surface reflections and 
diffraction . 

Assume that the interferometer optical system is alined. Then 
spiirioixs fringe shifts may be characterized by the angle <P which par- 
allel ray traces make with respect to the siirface of a model. The three 
basic ray-trace configurations with respect to a flat plate are illus- 
trated in figvire 4. The corresponding appearance of the surface of the 
plate for the three characteristic configurations utilizing first the 
test beam only, and then both beams, is shown in figures 5 and 6, re- 
spectively, for the condition of "no flow." Details of the observed 
fringe patterns are discussed in appendix C. In the presence of flow, 
as in a .wind tunnel, observed spurious interferences wo\xld severely com- 
plicate the interpretation of the Zehnder-Mach fringes, which are of 
primary interest. 

Although the complicated Interferences near the plate surface may 
be regarded as spurious because they tend to obscure the Zehnder-Mach 
fringes, nevertheless, they may serve useful pxirposes, namely, 

(1) Their existence provides a rapid and easy method for alining 
the light beam with the surface of the model. 

(2) The reflection fringes might be used to evaluate density dis- 
tributions in cooled boundary layers, that is, in situations where the 
light is refracted against the surface of the model. 

Methods for alining a light beam with a flat surface are described 
in references 3, 7, and 14. Alternatively, the alinement can be ac- 
complished by utilizing the spurious fringes, a simple magnifying lens, 
and a fine wire. As described in appendix C, the prominent surface 
phenomena are diagonal reflection fringes when <P < 0 (fig. 6(b)) and 
diagonal diffraction fringes when <p > 0 (fig. 6(d)). The two fringe 
patterns are inclined in opposite directions with respect to the sur- 
face of the plate. Alinement of the surface may be accomplished by al- 
lowing one end of the wire to touch the surface at the plate midspan 
plane. The image of the wire and the spurious interferences are ob- 
served by using the magnifier with the camera ground glass removed. The 
interferometer or the model is then rotated about the line formed by the 
intersection of the object plane and the model surface. Alinement of 
the plate is indicated by the "picket-fence" appearance of the spurious 
fringes at the surface (fig. 6(c)). Next, the ground glass is replaced. 
The fine wire may then be used for locating the ground glass at the de- 
sired image plane. Observing the point of wire contact with the plate 
surface also serves as a check upon the alinement procedvire. The entire 
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procedure is simple^ rapid^ and possibly somewhat more accirrate than 

the methods described in the references cited, although no data are 4 

available for comparing the methods. Residual misalinement resulting 

from application of the present method would definitely be well within 

the permissible misalinement limits set in reference 1. A similar allne- 

ment procedirre may be utilized for any surface which can be made parallel 

to the light beam. 


EVALUATION EQUATIONS FOR ONE-DIMENSIONAL DENSITY FIELDS, 

INCLUDING REFRACTION, END, AND CORNER CORRECTIONS 

In practice it is unlikely that any gas density field will be truly 
one dimensional. For example, in a wind tunnel the condition of one 
dimensionality will be violated in the vicinity of the spanwise ends of 
a model and by boundary layers formed on the wind-tunnel windows. Sim- 
ilar violations of other basic geometries are likely. With regard to 
interferometry, the optical effects introduced by boundary layers on the 
windows and at the spanwise ends of the model will be called end effects 
and corner effects, respectively. The following specific definitions 
will be adopted: 

(1) End effects - The optical effects caused by wind- tunnel-window 
boimdary layers which are involved in determining the free-stream 
density p 

00 

(2) Corner effects - The optical effects caused by boimdary layers 
in the vicinity of the spanwise ends of a model which are in- 
volved in determining the basic variable density field associated 
with the assumed geometry (p(y) in the present instance) 

The model may, or may not, be contained within a wind timnel. Thus, 
two cases of corner effects are of interest. If the model is contained 
within a wind tunnel, the model may, or may not, span the tunnel. If 
the model spans the tunnel, then the model and end- wall boundary layers 
will intersect. If the model does not span the tunnel, then the bound- 
ary layers may, or may not, intersect, depending upon how close the 
model is to the spanwise walls and other factors. In the present re- 
port the following two cases will be considered: 

(1) The model is contained within a wind tunnel and spans the 
tunnel. 

(2) The model is not contained within a wind tunnel. 

In many instances evalmtion eq.uations which apply to the other sittia- 
tions can be readily obtained by modifying the derivations presented 
herein . 
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Consider a basic one-dimensional density field p(y) . Assume that 
end and corner effects correspond to minor deviations from the basic 
assumed geometry. Major deviations would correspond to a different 
basic geometry. Thus^ it would then be appropriate to assume a dif- 
ferent basic geometry. An equation for computing the free-stream densi- 
ty Poo ^ wind tunnel and which Includes end corrections is derived 
in appendix D. Evaluation equations for p(y) are derived in appendixes 
E and F. These equations include corner corrections^ as well as re- 
fraction corrections. The previously listed case 1 is considered in 
appendix E, whereas case 2 is treated in appendix F. In these two cases 
evaluation equations are developed by extending the analysis reported 
in reference 1 to include corner effects. 


By means of Interferometry it is often possible to determine the 
free-stream density p^ within a wind tunnel from measurements per- 
formed outside the flow field. The appropriate experimental geometry 
consists of a wind tunnel of rectangular cross section bounded by win- 
dows at its spanwlse ends^ as shown in figure 7. Two adjacent spanwise 
reference holes are drilled through the spanwise walls or through a 
model which spans the tunnel. The hole denoted by r in figure 7 is 
connected to the extfernal atmosphere, whereas the hole denoted by w 
is connected to the timnel gas flow by means of a static tap. Appro- 
priate values of pressure, temperature, and density within each region 
are noted in figure 7. 

It is shown in appendix D that the free-stream density p is 
given by " 


where the window boundary layers are assvuned to be turbulent. The 
measurable fringe shifts and are associated with the 

density differences p. - p^, and p^ - p^., respectively, where Pj. is 

the density within hole r. The density p^ at a wind-tunnel, or 

model, wall which is parallel to the light path is also of interest 
and is given by 


Free -Stream Density p^ 




Further details concerning appropriate experimental and calculation 
procedures are given in appendix D. 


14 


NACA TN 3507 


Density Field p(y) (Model Spaxis Wind T\mnel) 

^ M /‘i Corner effects associated with a model contained within a wind 
f tunnel are considered in references 6 to 8. The appropriate geometry 
fi: consists of a wind tunnel of rectangular cross section hounded at its 

^ spanwise ends hy plane windows . A model whose shape is independent of 

the spanwise coordinate z completely spans the tunnel. The density 
field adjacent to the model surface is effectively a function of a sin- 
gle coordinate perpendicular to the surface, except at the spanwise 
ends, where the model-wall and window boundary layers intersect. A 
typical ray trace through the test section and relevant geometrical 
quantities are shown in figure 8. 

Evaluation equations corresponding to two particular assimed end- 
wall density distributions, namely, 

(1) An effective-average density 

(2) A power-law density distribution 

are derived in appendix E. Evaluation eqixations associated with distri- 
bution 1 are usefiil when the exact end distribution is unknown ♦ How- 
ever, in wind tunnels the boundary layers on the wind-tunnel windows are 
usually tiorbulent at the test section. Then, the end-wall density dis- 
tribution may be closely approximated by the power-law distribution 2. 

The evaluation equations corresponding to the two assumed end 
density distributions are as follows: 


Effective-average density: 

One-term approximation: = bQ 

p* = , 1+ H 



where 


Lg = L - I (Cl + ?2) 

^1 ^2 5 y ^0 


y H yo = yj) 
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Two-term approximation: M-r| = ^0 ^1^ 

7,2 

p* = 1 + N ^ l3? (4Lp - 6KL + " 


kP.Lg ■ 12 1 kp^L^ 


~ P2,vA/^ z ~ ^l^ 

^ p- A ^ /' 


where 

^ ^ “ 2 ^^1 '*’ ^2^ 

Lr = L - (^p + C2) 

^2 = ^ll|hiL2^ 


y = yo = yD - D 


D = ± |[(1 - 2K)L + ^2 - Ci]^>iLr 


^ = 


dJT 




where 

(2) 


= L 


Lg 

(2) _ 5z 


5 ( 2 ) 


£ 


s„ 


The designation of one- and two- term approximations refers to the 
number of terms of the assumed series expansion for refractive in- 
dex which are utilized (cf. ref. l) . The superscript (2) asso- 
ciated with Lg and indicates an approximation obtained by 

setting Cp = ^2 order to permit computation of the required 

quantities. Throughout the entire report the convention + or - is 
associated with p(y) increasing or decreasing, respectively. 
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Power- law distribution: 

One-teim approximation: = bg 


p* = 1 + N 


X 2 r- " ^2,w' 


% \ P. 

where 

Lr = L - (^1 + ^ 2) 





n + 1 


^1 


^1 “ ^2 " 8 ^0 

y 


y = yo = yp 


p = 


Two-term approximation: ” ^0 




n 


(5i + 5^) 


n + 1 


^1 = -r ^0 


Cg . ± i biL^ ^ 
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y = yo = - D 


D = ±|[(1 - 2K)L + ^2 - ?i]^iLr 


bl = 


dIT 

<3yD 




where 



L - 2^ 


( 2 ) 

1 


'’I by 


The approximations indicated by the superscripts in paxentheses 
are made for the same reason as in the preceding case. Calcula- 
tion procedures for both cases are similar to that given in ref- 
erence 1 . 


In order to include comer effects, additional quantities must be 
determined, namely, the window boundary- layer thickness 62 and the 


relative density- difference 



across the window boundary 


layer. Usvially the two quantities can be estimated or, if necessary, 
measured by some method other than interferometry. 


The present evaluation equations are more cumbersome than the 
corresponding equations presented in reference 1. Thus, in order to 
minimize computation, it appears advisable to calculate first the ex- 
treme magnitudes of the comer-effect contributions associated with 
the individual terms by comparing them with the corresponding terms 
in the equations given in reference 1. QirLte often, these Contribu- 
tions may be negligible compared with the primary and refraction ef- 
fects and, therefore, can be neglected, thus reducing the calculation 
to that indicated in reference 1 . 


Density Field p(y) (Model Not Contained within Wind Tirnnel) 

Corner effects associated with a model not contained within a wind 
tunnel are considered in references 4 and 5. The experimental geometry 
is similar to that for th^' preceding case except that the wind tunnel 
is absent. Again, the density field adjacent to the model is effec- 
tively a function of a single coordinate perpendicular to the surface, 
except at the spanwlse ends. A typical ray trace through the test sec- 
tion and relevant geometrical quantities are illustrated in figure 9 . 
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Evaluation eq.uatlons associated with an effective -average end 
density distribution are derived in appendix F. The resulting evalua- 
tion eq_uations are as follows; 

Effective-average density: 

One- term approximation: 


* 

P = 


PO 


= 1 + N 




where 

^ ^ ^ (^1 ^2) 

%, » L - 2(1^ - B^) 


H + 2(yp - By) 


y = yo = 

Two- term approximation: 
p* = 1 + N i bf (iLr + 3 


(1 - K)L - Lg + i ^ 2 ] 


where 

^ + 1 ^^1 ^2) 

= L - 2(1^ - 


Le = L 


^1 ^z + 2(y^ - Sp" 


^2 = ^z + 2Cyg - 6y) 
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yn = yo + s 

= * I 


CaJ 

OJ 

CD 


y = yo = ° 

D . * i[Lg - 2(KL + 6^ - JJJdiLr 


where 

4^) = L. 4^^ 




bi = 


dW 

dyo 




z(^> - I + 


‘z 2(yp - 6y) 

(2) 

^0 - 


V 


2tyD “ 


^ I 


Comments concerning the preceding case also apply in the present 
instance. The additional q.uantities which must be determined are 6y^ 

82^ and a. The model boundary- layer thickness 6y can be meas- 
ured directly from the interferogram. Often 82, and a can be 

measured from a supplementary interferogram obtained with the optical 
axis reoriented to be mutmlly perpendicular to the y-coordlnate axis 

and the original optical axis. The quantities ™ay be 

computed or may be meastired directly from the supplementary interfero- 
gram. Usually corresponds to atmospheric density. 


EXPERIMEOTAL APPLICATIONS 
Interferometer 

Three Zehnder-Mach interferometers possessing 6-inch-dlameter 
mirrors have been utilized at the NACA Lewis laboratory. The geometry 
of the latest design is shown in figure 10, and a photograph of the 
interferometer appears in figure 11. In a1 1 designs the mirror separa- 
tion ratio in perpendicular arms is approximately 2:1, thus permitting 
reorientation of fringes without defocusslng. Also, the commonly 
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practiced (refs. 15 to 17) but undesirable design characteristic^ if 
aberrations are considered, involving passage of the test beam through 
the final splitter plate is avoided. Optical specifications of the 
latest model follow: 

(1) Collimating lens: 30-Inch focal length, f/5 telescope 

objective 

( 2 ) Full-reflecting mirrors: 6-Inch diameter, 1-inch thick; 

aluminum-silicon monoxide coated; flat to O.IX sodium light 
(D line) 

( 3 ) Semireflecting mirrors: 6-Inch diameter, 3/4-inch thick; 

titanium dioxide coated; flat to O.IX sodium light (D line) 

( 4 ) Camera lens: Two of 20- inch focal length, f/4.5 Bausch and Lomt 

Tessars arranged face-to-face to foim an afocal system producing 
unity magnification 

The light-source unit provides for the use of incandescent, sodiiom, 
high-pressure mercury arc and magnesium spark light sources by the 
simple expedient of rotating a mirror inclined at 45° to the plane con- 
taining the light sources. The mercury arc lamp is eq.iu.pped for steady 
or flash operation. The magnesium spark unit and power supply are very 
similar to those described in reference 18. The spark unit produces 
about six times the luminous intensity of the mercury flash in the same 
exposure time (approx. 14(j.sec) . Sets of interference filters and 
apertures in sliding inserts are contained in the light- soiurce - 
collimator unit. The maximum permissible test-section span is 18 inches. 
A sliding insert in the reference path permits rapid insertion or re- 
moval of wind-tunnel-window compensator plates or other desired devices. 
The camera shutter and, if desired, a knife-edge mit for schlieren 
photography are located at the afocal point of the camera-lens system. 

All shutter controls are operated from outside the camera box by means 
of extension levers. The camera ground glass and film holder are part 
of a standard 5-by-7 plate-back camera and are adjusted for focussing 
by means of the camera bellows. A 35-milllmeter film-adapter unit is 
also available . 

The afocal camera- lens system was selected to produce unity magnifi- 
cation of all object planes throughout the test-section span. Thus, no 
correction for image magnification is required. The actual measured 
magnification as a function of the object plane is listed in the fol- 
lowing table, where the entire available test-section span is 18 inches: 


SpazLwise 

Magnifi- 

distance^ 

cation 

in. 


4.5 

0.9984 

9.5 

.9979 

14.5 

.9985 


3378 


3378 


NACA TN 3507 


21 


The interferometer light source and camera subassemblies are quickly 
removable and replaceable. The interferometer is supported by adjustable- 
height legs eqiiipped with castered wheels; trunnions permit rotation of 
the Instrument when supported by a yoke. The aforementioned accessories 
permit movement^ installation, and adjustment of the Instnment at a new 
facility a ll in one day. Moreover, a complete set of boundary-layer 
data (24 interf erograms , foiar flow conditions) has been obtained in not 
much more time than that required to obtain the desired wind-tunnel 
conditions (about 2 hr). 

A discussion of three typical experimental applications of gas-flow 
interferometry follows. Novel exp lanations for dlscrepsuicles between 
theory and experiment are proposed. 

The applications are reported in the sequence in which the experi- 
ments were performed in order that successive improvements in technique 
might be noted. 


Supersonic Air Flow Along a Flat Plate 

Laminar -boundary- layer density profiles associated with supersonic 
air flow along an insulated flat plate and reported herein were calcu- 
lated from data utilized in refer ence 10. The present density calcula- 
tions differ from those reported in reference 10 in that the evaluation 
equations derived in reference 1 and in the present report are applied. 
Appropriate data for the present calculations were 

L = 3.60 in. 

X = 2.15x10“^ in. (mercury green line) 

K = 1/2 

k = 0.117 cu ft/slug 
M « 2.04 

. «o 

0.4X10® < Re <1.0X10® 

X 


X =! 2.5 in. 

where is the free-stream Mach number, and Rex is the Reynolds 

number based on the distance x downstream of the leading edge of the 
plate . 
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Following estimation of corner effects,, density profiles were com- 
puted according to the two-term equation presented in reference 1. The 
computed profiles are shown in figinre 12(a)- as a function of the dimen- 
sionless distance ^V'Re^. The corresponding theoretical density pro- 
file, according to Chapman and Rubesin (ref. 19), is indicated for com- 
parison with the experimental data. Velocity profiles u* as a function 

of ^VReZ, where u* = were calculated from the eqiiation 

X * 


u* = 


Xr - i)K 




and are shown in figure 12(b) for comparison with the theoretical pro- 
file predicted in reference 19. The theoretical total-temperature ratio 

T* = indicated in reference 19 was assimied for computing u* 

from ejqperimental data. Values of the section skin-friction coefficient 
C were calculated by numerical integration of the modified Karman 

momentum equation 



The computed coefficient C^ as a function of Re^ is shown in figure 
12(c) for comparison with theoretically predicted results for laminar 
flow. 


Density profiles were also computed according to the one- and three- 
term approximations derived in reference 1 and the present two-term ap- 
proximation, which includes comer effects. The resulting density pro- 
files for the largest Reynolds number are shown in figure 13 in order to 
illustrate the comparative importance of refraction and corner effects. 

The recomputed density profiles shown in figure 12(a) exhibit a 
greater tendency to collapse into a single profile than do those shown 
in reference 10. Also, there is a general shift of the profiles so that 
the value of p* associated with a given value of y is larger than 
was the case in reference 10. The experimental density and velocity 
profiles are very similar to those exhibited in reference 8. A system- 
atic disagreement between theory and experiment is evident. The recom- 
puted values of are in better agreement with theory than the re- 

sults presented in reference 10. The difference between experimental 
and theoretical values of varies from 0 to 30 percent, with the 

greatest relative difference associated with the largest Reynolds nmber. 
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The maximum refraction correction varied from 4 percent (at 
“ 0.483X10®) to 10 percent (at = 0.989X10®) of the relative 

TTiayi TTiuTn density difference (p^ - P )/p^* The following additional 

items are important in detennlning the validity and error of the com- 
puted density profiles (cf. ref. l) : 

(1) Appllcahility criterion: Determines when the two-term approxi- 

mation equations are valid 

(2) Apparent-ray- trace-crossing criterion: Specifies the condition 

for which the interferogram is free of superimposed imagery 


(3) Light-source criterion: Specifies magnitude of mlsalinement 

or lateral extension of light source for which the assumption of 
an axial point source^ hence the proposed evaluation equations, 
remains valid 

(4) Series remainders: The errors in p* and D resulting from 

termination of the series expansion for refractive index 

(5) Comer effects 

(6) Error in measuring fringe shifts and free-stream density p 

OD 

The computed extreme value associated with each of the preceding 
items was as follows: 


(1) Applicability criterion: 


-2X10^ in. <d%/dy2 < 2X10^ in.-^ 

-2X10^ in."^ < d^N/dyp < 0 

(2) Apparent-ray- trace-crossing criterion: 

|dD/dy| <1.0 (criterion) 
|dD/dy| <0.1 

(3) Light-source criterion: 

(q/f)2« 1.0 


(criterion) 

(experiment) 


(experiment) 


(q/f)2 < 4X10"^ 

(p/f)2 « 1.0 

(p/f)2 < 10"^ 


(criterion) 

( experiment ) 
( criterion) 
(experiment) 
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(Quantities p and q. are lateral coordinates of points in the plane 
of the light source . Coordinate p corresponds to and cor- 
responds to y. The origin (p^q) = (O^O) lies on the optical^ or , 
axis . ) 

(4) Series remainders: Corresponding to the two-term approximation^ 

[RgCP*)! ^0.002 

|R2(D) I = 2X10"^ in. 

(5) Comer effects: The errors introduced by neglecting corner 

effects were found to be 


Ap* < 0.01 
AD < 0.0001 in. 

(6) Error in measuring N and p^: For 'the interferograms pre- 

sented in reference 10, dj^ * 0.01 inch. The fringe shifts were care- 
fully measured. Therefore, by assuming that the data in figure 2(b) 
apply and by utilizing the error emresslon for orientation parallel to 
grad p, it is estimated that |ANy^| - 0.02 in the region (adjacent to 
free stream) associated with the greatest disagreement between theory 
and experiment. In 'the indicated region, N « -0.5. Thiis, the standard 

deviation of p* caused by the error in measuring N was /Ip* * 0.001. 
The error in measuring p^ is unknown. 

All criteria were satisfied. In fact, it is quite unlikely that 
the liglit-source criterion would ever be violated unintentionally. The 
values given for the series remainders are much less than the values 
originally computed according to the remainder formulas given in ref- 
erence 20. The corner effect and random error in p* caused by the 
error in measuring N were negligible. 

The preceding errors are insufficient to account for the systematic 
disagreement between theory and experiment, especially in the region 
adjacent to the free stream where refraction corrections are negligible. 
Remaining sources of possible significant errors are 

(1) Error in measirring free-stream density p^ 

(2) Flow transition 

(3) Transverse contamination 

(4) Incorrect assumptions regarding the boxmdary- layer geometry 

(5) Nonuniform heating of the wind-tunnel windows 
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(1) The disagreement between theory and experiment indicated in 

figure 12(a) is systematic. The error in measuring is likely to he 

random. A systematic error in corresponds to flapping the profile 

about the ordinate value p* = 1 with the profile at y = y^ held 

fixed and the end of the profile at y = 0 free. Agreement between 
theory and experiment is not attainable by this procedirre. 

(2) The disagreement with theory is not confined to the present 
results. A similar disagreement is apparent in the results of inter- 
ference measurements presented in references 3 and 8. Although the 
measured density profiles exhibit certain ctiaracteristics associated 
with transition, in that they fall between theoretical laminar and 
turbulent profiles, the profiles presented in references 3 and 8 tend 
toward agreement with laminar theory for increasing Keynolds number 
associated with increasing x. This is contrary to expectations if 
transition were a factor. Moreover, in reference 8 shadowgraph photo- 
graphs were used to confirm laminarity. Thus, flow transition is ex- 
cluded as a cause of the systematic disagreement. 

(3) Transverse contamination (refs. 21 and 8) in one sense, at least, 
does not appear to be the source of the error because the disagreement 

is similar both for transient (ref. 8) and steady-state (present re- 
port) experiments. For transient conditions the temperature in the 
corners is lower than that of the plate, whereas for steady-state con- 
ditions the corner temperature should be higher. Thus, the resulting 
density profiles should be somewhat different in the two cases. A 
similar conclusion is indicated in reference 8. 

(4) Although the present corner-effect analysis is admittedly crude, 
it should permit calcxilation of the effect to, at least, the correct 
order of magnitude. The interaction of the plate and window boundary 
layers is not abrupt at right angles, as was assumed in the analysis; 
but, rather, a smooth fairing together of the boundary layers may be 
expected. However, if the boundary- layer cross-sectional shape is as 
shown in figure 14(a) or reference 22, then the resulting correction is 
opposite to that which would improve agreement between theory and ex- 
periment. If the cross-sectional geometry is as shown in figure 14(b), 
then the disagreement can be completely accounted for. Prevalent flat- 
plate laminar -boundary- layer theories (e.g., ref. 19) are based upon 

the assumption that the plate span is infinite. However, transverse 
bulging of boimdary layers on the walls of supersonic wind runnels has 
been found (ref. 23), and it appears not \mlikely that similar effects 
may occ\ir in the boundary layer on an interposed flat plate. Ulti- 
mately, the bulge, which apparently arises from transverse pressure 
gradients, might be expected to decay in the constant-area section 
downstream of the leading edge of the plate. Agreement between theory 
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and experiment tends to improve with increasing x (refs. 3 and 8) in 
concordance with this argument. The present results indicate improved 
agreement between theory and experiment with decreasing free- stream 
density at a fixed value of x^ which corresponds to decreased bulging 
with decreasing density at fixed x. The effect is confirmed by probe 
measurements (ref. lO) at three spanwise positions. 

(5) Finally^ temperature variations throughout the wind-tunnel win- 
dows are accompanied by variations of window thickness and refractive in- 
dex. Temperature variations may exist at the window surface adjacent to 
the flow^ thus causing nonuniform heating of the windows. For exan^le^ 
measured surface -temperature differences of 2 ^ C between adjacent re- 
gions of laminar and turbulent flow on a flat plate are reported in 
reference 8. Section surface-temperature differences on wind-tunnel 
windows are undoubtedly considerably less than 2 ^ C because the boundaiy 
layer is then entirely turbulent. However^ from a simplified analysis 
presented in appendix it is found that a transverse temperature var- 
iation of only 0.1^ C would introduce an absolute error in p* amount- 
ing to 0.005^ which is on the verge of being appreciable. Further in- 
vestigation of the magnitudes of such temperatmre variations is necessary 
before any definite conclusions regarding their importance can be 
reached. In any events the resulting error can be minimized by choos- 
ing window glass having the minimum possible effective combination of 
coefficient of thermal expansion and temperature coefficient of refrac- 
tive index and by keeping the ratio of window thickness to wind-tunnel 
span as small as is mechanically feasible. If necessary^ a crude cor- 
rection for the effect can be performed by obtaining an interferogram 
immediately following termination of flow. In experiments not reported 
herein the effect of nonuniform window heating has appeared as a false 
boundary layer following termination of flow. 


Subsonic Air Flow Along a Thick Flat Plate with Blunt Leading Edge 

Boundary- layer density profiles associated with subsonic air flow 
along a thick flat plate with blunt leading edge were determined by 
means of interferometry for comparison with similar profiles obtained 
by two other methods of measurement^ namely^ X-rays and a total-pressure 
probe (ref. 24). A sketch of the model is shown in figure 15. A com- 
plete description of the wind tunnel and model is given in reference 
24. Data pertinent to the interference calculations were as follows: 

L = 3.806 in. 

X = 2.15x10"^ in. 

K = 1/2 
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k = 0.117 cu ft/slug 
0.55 < M < 0.78 

«o 

10^ < Re^ < 5X10^ 


X = 2.74 in. 

Six interferograms were obtained for each of four Reynolds niomber - 
Mach number combinations. Typical Interferograms obtained by using im- 

f iltered and filtered mercury light ( X. = 5461 A) are shown in figures 
16(a) and (b)^ respectively. Averaged density profiles for each of the 
four conditions are shown in figiire 17 as functions of the dimensionless 
distance y/0^ where the momentum thicknesses 9 were computed by 
numerical integration of the Karman momentum equation. No attempt has 
been made to compare the experimental results with theory. However^ for 
at least one flow condition (Rex = 0.443X10^j M, = 0.55) all six inter- 
ferograms indicated densities in the boundary layer adjacent to the 
free stream which were slightly greater than p . This might arise from 
the fact that the air flow made a 90° tiom in proceeding from the plenum 
chamber to the test section. That such a condition may be possible is 
shown in reference 25. 

In the present experiment the interferometer was used for deter- 
mining the free-stream density and the density p^ at the surface 

of the model, as well as the density field p(y). Because the plate 
spanned the wind t\mnel, it was possible to drill two l/4-inch-diameter 
reference holes spanwise through the plate at the chordwise location 
(fig. 15) at which boundary- layer measurements were desired. Cross 
hairs were located within the holes at the wind-tunnel midspan plane, 
which was the selected object plane. The wall static hole w was con- 
nected to the tunnel air flow by means of a 0.005-inch-diameter hole 
which was drilled perpendicular to the plate surface near the midspan 
plane. The reference hole r was connected to the external atmosphere 
by means of a small hole which was drilled thro\;igh one of the wind- 
tunnel windows. The procedure outlined earlier in the present report 
and detailed in appendix D was used for determining p and p^. 

The preceding technique possesses several desirable attributes, 
namely: 

(1) The flow is completely undisturbed by measuring Instruments. 

(2) All aerodynamic pressure-measuring Instrumentation can be 
eliminated. 
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(3) Additional req.uired instrumentation consists only of a barom- 
eter, a thermometer, and a thermocouple for measuring the model 
temperature • 

(4) The entire density field, including free-stream and model wall 
densities, is recorded instantaneously on a single photographic 
negative • 

(5) Instantaneous density profiles based on simultaneous values 
of p are obtained. 

00 

(e) Because is determined from the interference patterns ap- 

pearing in the reference holes , the problem of spurious fringe shifts 
caused by the model surface is circumvented. 

(7) Temperature gradients within the model are indicated by the 
interference patterns appearing in the reference holes. Suitable modi- 
fications of the basic technlq.ue can usually be made in order to ac- 
count for such gradients. 

Three different methods of determining the average wall-to-stream 

density ratio p* = p /p were compared. The ratio was determined 
w «> 

from the Interference data by means of the technique already described, 
from pressure measurements by using the equation 



1 

1 + 0.2tiM^ 


where the value of the temperature recovery factor t] was assumed to 
be 0.88, and from the pressure and temperature measurements by using 
the equation 


_ ('^t^oo 1 

Pw ~ o 

Tw 1 + 0.2M^ 


where (T^)„ was measured in the plenum chamber. 

The averaged density ratios corresponding to each of the flow con^ 
ditions are compared in the following table; 
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Condi- 

tion 

Rex 

M 

m 

p* from 

pressure 

probe 

p* from 

pressure and 

temperature 

probes 

p* from 

inter- 

ferometer 

1 

0.179X10® 

0.72±0.04 

0.916*0.011 

0.923 

0.915*0.016 

2 

.442 

.7831.008 

.902*. 006 

.905 

.9011.007 

3 

.443 

. 5501.008 

.9491.005 

.953 

.9461.003 

4 

.952 

.782*. 008 

. 9021.006 

.906 

• 905± • 002 


The maximum disagreement between the values of is 0*9 percent. 

Corresponding values of p* obtained by the three methods agree to 
within the standard deviations of p*. For the pressure measurements^ 
the pressure fluctuations were of the order of +0-05 inch mercury for 
conditions 2 and 3 and ±0.1 inch mercury for conditions 1 and 4. It was 
apparent from the observed interference patterns and the density pro- 
files computed from them that variations of p* resulted principally 
from flow fluctuations rather than inadequacies of the three methods 
for steady-state measurements. 

The fringe pitch d^ was necessarily varied from 0.235 millimeter 

(condition l) to 1.45 millimeters (condition 4) in order to obtain the 
zero-order interference fringe in both holes simultaneously (cf . appen- 
dix D) . The corresponding error AN^ varied from about 0.01 to 0.005^ 
respectively. The values of are considerably less than that of 

AN^ namely^ AN » 0.04 for all flow conditions^ where AN represents 
the standard deviation associated with all interferograms corresponding 
to a given flow condition. The latter value includes the effect of the 
error (±0.0002 in.) in locating the model surface with respect to the 
interference pattern. Because AN was from four to eight times great- 
er than AN^^ it is concluded that the scatter of the fringe-shift pro- 
files must have been caused primarily by mass-flow fluctuations. In the 
present instance only steady-state conditions were of interest. In the 
presence of the minor flow instabilities^ steady-state conditions were 
approximated with sufficient accuracy by averaging six interferograms. 

Extreme vaJLues associated with the various criteria and errors 
were as follows: 

(l) Applicability criterion: 

-2X10^ in.“^ < d%/dy^ < 2X10^ in.'^ (criterion) 

-0.5X10^ in. < d%/dyp <0 


(experiment) 
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(2) Apparent-ray-trace-crossing criterion: 

|dD/dy| <1.0 (criterion) 

IdO/dyl <0.007 (experiment) 

(3) Series remainders: For the two- term approximation, 

|k2(p*)| < 10-4 
|R2(D)| < 0.0001 in. 

(4) Corner effects: The maximum errors introduced by neglecting 

corner effects were 


Ap* < 0.002 

AD « 0 

All criteria,; including the light-source criterion,, were satisfied 
and the preceding errors were negligible. The maximum refraction cor- 
rection varied from 1 percent (at Re^ = 0.179X10^) to 3 percent (at 

Re^ = 0.952x10^) of (p^ - P^)/p^- Side-wall boundary-layer density pro- 
files were determined by means of a total-pressure probe and were foiind 
to satisfy a l/9 power-law distribution. ThuS; the evaluation equa- 
tions associated with a power-law end-density distribution (n = 9) were 
used for computing the corner effects. 


Free Convection of Air Around a Heated Horizontal Cylinder 

Free convection of air around a horizontal circular cylinder has 
been investigated theoretically (ref. 26). Experimental measurements 
(refs. 26 and 5) have given good agreement with theory. Thus, the tem- 
perature profile beneath a cylinder can be regarded as known and can, 
therefore, be used for testing the analysis herein. The cylinder problem 

(1) Serves as a test of the evaluation equations for a model not 
contained within a wind tunnel 

(2) Pemiits comparison of results obtained by treating the temper- 
ature distribution as a one-dimensional distribution and as an 
axially symmetric distribution 

(3) Permits comparison of the present experimental results with 
previous theoretical and experimental results 
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A Duralumin cylinder 3 inches in diameter and 10 inches long was 
utilized in the experiments. Further details of the cylinder construc- 
tion and the experimental procedure are discussed in appendix H. 

Sets of interferograms were obtained corresponding to the three 
optical orientations depicted in figirre 18^ namely^ 

(1) Light propagation parallel to the cylinder axis (one- 
dimensional distribution with corner effects) 

(2) Light propagation perpendicular to the cylinder axis at the 
midlength plane (axially symmetric distribution) 

(3) Light propagation perpendicular to the cylinder axis at an end 
plane (depicts corner distribution for computing corner 
corrections) 

Typical interferograms characteristic of the three optical orientations 
are shown in figure 19. Orientation 1, shown in figvire 19(a), permits 
evaluation of the vertical temperature distribution beneath the cylinder 
according to the procedure for a one- dimensional field associated with 
a model not contained within a wind tunnel. Orientation 2, shown in 
figiare 19(b), permits evalxoation of the same field according to the 
axially symmetric method. Orientation 3, shown in figure 19(c), pro- 
vides Information concerning the density distribution at the ends of 
the model. 

Interferograms were obtained for several wall-to-ambient tempera- 
ture differences © in the interval 24° F < 0 < 450° F by using orien- 
tation 1 with K = 1/2 and l/3 and, also, by using orientation 2. Ac- 
cording to the theory of Hermann (ref. 26) all temperature profiles 
wiiich satisfy the conditions 


where Gr is the Grashof number, 5 is the boundary- layer thickness, 
and r is the cylinder radius, should collapse into a single profile 
for any given azimuth angle when the dimensionless temperature ratio 


and g(x) is an azimuth function which possesses the value 0.760 beneath 
the cylinder. In the present experiments the quantities Gr and &/r 
were foimd to satisfy the Inequalities 


10^ < Gr < 10^ 
6/r « 1 



7.4X10^ < Gr < 2. 3X10 
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when the kinematic viscosity v was based on the cylinder wall temper- 
ature, and 

0.25 < 5/r < 0.30 

Therefore, except for the fact that S/r was comparatively large, the. 
conditions of the experiment satisfied the req.uirements of the theory. 

Dimensionless temperat\are profiles obtained for orientation 1 with 
K = 1/2 and computed according to the two-term approximation with • 
corner corrections are compared in figure 20(a) with the theoretical- 
profile predicted in reference 26. (incidentally, the dimensionless 
coordinates of the theoretical profile are not given in ref. 26 but can. 
be obtained from table I of ref. 27.) Profiles obtained for K = 1/3 
and computed, first, according to the Wachtell-DeFrate method (refs. 2.^ 
and 3), which neglects corner effects, and then according to the present 
two-term approximation are compaxed with the theoretical profile in 
figures 20(b) and (c), respectively. Profiles obtained by using orien- 
tation 2 and computed according to Wachtell's simplification (ref. 28) 
of Von Voorhis' evaluation procedure (ref. 29) for axially symmetric 
distributions are compared in figure 20(d) with the theoretical profile. 


The kinematic viscosity 
, 3 . 


V, which is reqiiired for computing the 


Grashof number, Gr = ^ (I^ - l), where g is the gravitational accel- 

8v^ 

eration and T^ is the wall-to-ambient temperatiire ratio, was in all 

The effect of basing 


cases based on the cylinder wall temperature T.^.. 
V on the ambient temperatiore T^, rather than 
ure 21. 


T , is shown in fig- 


Agreement between the experimental and theoretical temperature pro- 
files is generally very good and somewhat better than that indicated in 
reference 5. This should be expected because, in the experiments re- 
ported in reference 5, S/r was somewhat greater (s/r » 0.5) than in the 
present experiment and refraction was apparently neglected. In refer- 
ence 5, the kinematic viscosity v was based on (T^ + T^)/2, which was 

assumed to be constant for all values of 0 . By comparing the tempera- 
ture profiles in figures 20(a) and 21 with those in reference 5^ it is 
apparent that basing v on (T^ + 'the cause of the dis- 

crepancy between theory and experiment which was obtained in reference 5. 

The maximum refraction correction varied from 0.2 percent (for 
0 = 24^) to 2 percent (for 0 = 450^) of (p^ - correspond- 

ing maximum corner corrections were found to vary from 0.8 percent (for 
© = 24°) to 8 percent (for 0 = 450^). Therefore^ refraction and corner 
corrections were included in computing the temperature profiles. 
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LO 


Extreme values associated with the various criteria and remainders 
were as follows: 


(criterion) 

(experiment) 

(criterion) 

(experiment) 


(1) Applicability criterion: 

-0.8X10^ in. “2 < d%/dy§ <0.8X10^ in."^ 

-1.2X10^ in. -2 < d2N/dy2 <0 

for K = l/2 and 0= 450.1° F; 

-4X10^ in. "2 < d%/dy§ < 7X10^ in."^ 

-O.lxlO^ in.”2< d%/dy^ <0 

for K = 1/3 and © = 202.0° F. 

(2) Apparent-ray-trace- crossing criterion: 

|dD/dy|<1.0 (criterion) 

|dD/dy|<0.006 

for K = 1/2 and 0 = 450.1° F; 

I dO/dy |< 1.0 

I dO/dy I < 0.0006 

for K = 1/3 and 0 = 202.0° F. 

(3) Series remainders: For the two- term approximation, 

R2(p*) < 0.008 
R 2(D) t 0.003 in. 

R2(p*) = 0.8X10"^ 

R2(D) = 0.8X10"^ in. 


(experiment) 

(criterion) 

(experiment) 


for K = 1/2 




). for K = 1/3 


The light-source criterion was easily satisfied. The applicability 
criterion was violated at the largest values of 0 for K = 1/2. The 
criterion, therefore, must be conservative, as was indicated in refer- 
ence 1, because the corresponding temperature profiles in flgvires 20(a) 
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and (c) are all in good agreement with theory. The distortion remainder 
of 0.003 inch for K = l/2 is only about 0.7 percent of the boundary- 
layer thickness and^ of course^ represents the maximum possible value 
attainable by R2(D) . 


CONCLUDING REMARKS 

In the three experiments, the ratio of wind- tunnel -wall boundary- 
layer thickness to wind-tunnel span possessed the following maximum 
values : 

Supersonic flow: 262/L » 0.06 

Subsonic flow: » 0.16 

Free convection: 252/L « 0.025 

Corner effects were negligible in the supersonic- and subsonic-flow 
experiments but were quite large in the free -convection experiment, 
especially for large values of 0 . From these results, it is concluded 
that relatively thick wlnd-tiannel-wall boundary layers are likely to 
introduce negligible or small corner effects (although end effects are 
usually significant), whereas relatively thin end boundary layers are 
likely to introduce relatively large corner effects when the model is 
not contained in a wind tunnel. Therefore, the existence or nonexist- 
ence of windows bounding the model is far more Important in determining 
the significance of corner effects than is the relative end-wall- 
boundary-layer thickness. 

The applicability and especially the ray-crossing and light-source 
criteria were effectively satisfied in all three experiments. Only the 
distortion remainder in the free-convection experiment appeared to be a 
possible significant residual error among the systematic errors which 
were evaluated. 

Knowledge of the random error AN^ as a function of dj^ (fig. 2(b)) 

should prove useful in determining the error in measuring time fluctua- 
tions of N in tiirbulence studies (refs. 6 and 30) because AN^ is 
exactly that error which would be of interest. 

In reference 31 it is shown that the usefulness of interferometry 
in studies of cooled boundary layers is theoretically limited by re- 
fraction of light against the model surface. However, if the model sur- 
face is a reflecting surface, then light is reflected from the surface 
rather than absorbed by the surface as is assumed in reference 31. A 
typical interferogram shown in figure 22 and obtained by using the 
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hot-plate model described in reference 1 illustrates the reflection 
effect. At the cool s\arface (the lower sirrface in the interferogram) 
the interference fringes do not terminate but^ rather^ are continuous^ 
appearing as reflection fringes below the surface. The surface is in- 
dicated^ approximately, by the horizontal line passing through the 
vertices of the "parabolic" portion of the fringes. A method for 
svS'luatlng the temperature distribution adjacent to the cool surface 
(in the example presented) might be developed by extending the analysis 
presented in reference 1 to include reflection and the phase shift 
accompanying reflection (ref. 32) . Success of the method would 
obviate the difficulties described in reference 31. The method might 
also be applicable when the light is refracted away from the surface 
by employing an off-axis light source expressly for the purpose of ob- 
taining reflection fringes. Consideration of off-axis light sources 
is included in the analysis contained in reference 1. 

Although the interference method permits a vast amount of data to 
be obtained in a very short time, conversion of the data contained in 
the interferogram into numerical data is at present a time-consuming 
process. The usefulness of the Interference method coiild probably be 
increased most significantly by the development of an automatic, or 
semiautomatic, instrument for converting the data contained in the in- 
terference pattern into numerical data. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 13, 1955 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

a perpendicular distance from outer hyperbolic- shaped extremity of 

corner boundary layer to intersection of hyperbola asymptotes 

(fig. 9) 


to 

value of refractive index associated 
for any given ray trace 

with 

T) = 0 

and 

o 

II 

tl 

refractive-index gradient associated 
for any given ray trace 

with 

T] = 0 

and 

II 


b<y refractive-index coefficient associated with t] = 0 and y = yQ 

for any given ray trace 

Cx section skin-friction coefficient 

c^ light-path coefficients (ref. l) 

c* light-path coefficients (ref. l) 

fringe-shift coefficients (ref. l) 

D optical distortion 

d interference -fringe pitch 

Gr Grashof number 

H final refraction displacement perpendicular to z-axis and within 

density field p(y) 

Zr 

K fraction of test- section span^ K = 1 - — 

JL 

k specific refractivity 

L model span 

Lg effective model span in wind tunnel, Lg = L- - 2 (^]_ + ^ 2 ) 

T.^ ^) approximation of ~ ^ 


3378 


R/.P<7 


NACA TN 3507 


37 



effective span of model not contained in wind tunnel^ 

^ ^ f (^1 ^ 2 ) 

l(2) 

approximation of 

S: Le 

o: 

effective wind-t\mnel span in determining end correction, 
^ ■ n + 1 

Le 

effective span of model not contained in wind tunnel, 
Lg = L - Jzi 


reduced model span in wind tiumel, Lpj = L- (Ci + ^g) 

4^> 

( 2 ) ( 2 ) 

approximation of Lp,L^'=L-2^£' 

r* 

reduced span of model not contained in wind tunnel, 
Lr = L - 2(1^ - 


extent of corner boundary layer in z -direction 

M 

Mach number 

m 

lateral magnification 

N 

interference-fringe shift in units of fringes 

% 

observed interference-fringe-shift profile^ Np = Nj)(yj)) 

Nw,r 

interference-fringe shift associated with density difference 
Pw “ Pp 

N „ 

interference- fringe shift associated with density difference 
- Pr 

n 

integer 

n 

exponent of power-law density distribution 

P 

absolute pressure 

R 

gas constant 


remainder associated with power- series expansion of \x to two 

terms 
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Rex 

s 

T 

Tt 

t 

u 

u* 

X 

y 

yo 

y(2) 

z 

2R 

P 

P 

r 

A 

A 

5 


Reynolds mmber based on length x 

physical length 

static temperature 

total temperature 

window thickness 

velocity 

velocity ratio, u* = u/u^ 

coordinate perpendicular to optical axis 

value of X associated with ^ = z = 0 for any given ray trace 

coordinate perpendic\lLar to optical axis 
y- coordinate value of measured fringe shift 

y- coordinate value of ray trace when t| = H 
y- coordinate value of ray trace when tj = 0 
approximation of y^, y^^) = y^^ 
coordinate parallel to optical axis 

distance from origin of z to selected object plane 

projected angle formed by ray trace and optical axis in yz-plane 
temperatirre coefficient of refractive index 

ratio of specific heat at constant pressure to specific heat at 
constant volume 

error 

increment 

boimdary- layer thickness 

boundary-layer thickness associated with density distribution 
P = p(y) 
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■boundaxy- layer thickness associated with refractive index distri- 
bution (ig 


ray-trace coordinate parallel to optical axis, ^ = z 

assumed z component of ray trace in corner adjacent to initial 
boundary z = 0 (fig. 8) (model contained in wind tunnel) 


approximation of 



assixmed z component of ray trace in corner adjacent to initial 
boundary z = 0 (fig. 9) (model not contained in wind tunnel), 

_ a^ 

^1 = ^z + 2(y^ - 5y) 

assumed z component of ray trace in comer adjacent to final 
boundary z = L (fig. 8) (model contained in wind tunnel) 

assumed z component of ray trace in comer adjacent to final 
boundary z = L (fig. 9) (model not contained in wind tunnel), 

7 , a" 

^2 - + 2(y^ - 5^) 

ray- trace coordinate parallel to y-coordinate, tj =|y - yQ| 
absolute temperature difference, ® - T 


absolute temperature difference, 6 = T - T. 


wave length of light 
refractive index 

refractive -index distribution adjacent to model, p = p(y) 
refractive index of wind-t\xnnel windows 

refractive-index distribution adjacent to model as function of 
coordinate t] 

value of refractive index associated with tj = 0 and y = yQ 
for any given ray trace 
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^^2 

1^3 

V 

V 

C 

P 

PQ 


w 


refractive -index distribution adjacent to wind-tunnel windows, 

H 2 = 1^2(2) 

refractive-index distribution in vicinity of corners of models 
d3 = PsCy^z) 

kinematic viscosity 

integer 

ray- trace coordinate parallel to x- coordinate^ ? sjx - Xq| 
density 

value of density associated with t) = 0 and y = y^ for any 


given ray trace 
density ratio ^ P* = Po/P« 

value of density ratio at model surface^ p* = P^/p^ 
a integer 

cp angle formed by ray trace and surface of flat plate 

\|r optical-path length 

Subscripts : 
a atmosphere 

integer 


n 


w 

V 

a 

2^w 


reference value within hole spanning wind tunnel and connected 
to external atmosphere 

model surface or wind-tunnel wall parallel to model surface 

integer 

integer 

wind-tunnel window surface 
ambient^ or free-stream^ value 
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Superscripts: 

mean value 

' refers to image space associated with camera lens 


CD 


g 
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APPENDIX B 


RANDOM ERRORS IN FRINGE-SKEFT MEASUREMENTS 


Random errors in measuring N were determined from interf erograms 

obtained with monochromatic light ( X = 5461 A) . Two observers performed 
the measiirements . The pitch of two adjacent interference maximums or 
minimums was measured five or 10 times ^ as noted in figure 2. Scatter 
of the resiilts in figure 2 proved to be independent of the observer and 
independent of whether the measurements involved Interference maximums 
or minimums. However, contrast variations influenced the scatter, the 
computed values of Ad^ being appreciably less for high-contrast inter- 

f erograms than for low- contrast interf erograms . The effect of the ratio 
of widths of interference maximums to minimimis was not explicitly in- 
vestigated, although the ratio varied from about 0.4 to 0.8. 

It should be noted that the experimental results are based upon 
measurements made perpendicular to straight parallel fringes. When p 
is an arbitrary function of y, the fringe- shift measurements are made 
diagonal to cvi^ed fringes. Because the effective pitch in the x'- 
direction is the same when p = or P = p(y)^ the value of Ad' 

associated with measurements between the inclined fringes may be some- 
what greater than the values presently reported because of the diffi- 
culty of exactly locating maximums and minimums along a diagonal path 
of approach. In addition, if the second derivative of p is large, 
then interference maximums and minimums will be displaced from the 
centerline of the bright and dark fringes, respectively, thus making 
settings on mayiTnums and minimums indefinite. However, the fact that 
the normal pitch of the Inclined bright and dark fringes is reduced may 
provide a compensating factor which tends to reduce Ad'. 

The two-term approximation = b^ + bpT)) equations for p* and 
D are, neglecting corner effects. 



D = ± I (1 - 2K)b^L2 


where 


_ ^ ^ 

1 dyj) L 
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The preceding equations are derived in reference 1. The corresponding 
primary maximum random errors of p* and D are given by 


Ap* = (1 - p*) 



AD ^ m; . , al , ax 

D ~ dlj/dyp (1 - 2K) L X 


The importance of the error AN/n is of interest. Thus^ as p* 1, 

1 - p* 0, and AN/n -*■ Therefore^ AN/n is the primary error. 
However, as p* 0, 1 - p* -»■ 1, and AN/N -*■ 0. Then Z^/n becomes 

less significant and^ in particular^ may be insignificant compared with 
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APPENDIX C 


SURFACE INTERFERENCE EFFECTS 

The primaiy Zehnder-Mach Interference pattern resvilts from the 
superposition of the test heam^ which traverses the test section con- 
taining the model, and the reference beam, which circumvents the test 
section. The observable forms of interference induced by a flat plate 
contained within the test section and oriented nominally parallel to 
the light beam appear to be primarily comprised as follows: 

(1) Test beam plus reference beam (Zehnder-Mach interference) 

(2) Diffraction by plate edge nearest light source 

(3) Reference beam plus reflected beam 

(4) Reference beam plus light diffracted by plate edge nearest 
image plane 

Other so\rrces of interference are not vital to the present discussion. 

The first two mentioned forms of interference are quite common. The 
third form is similar to Lloyd's mirror interference. No mention of 
the foiirth form has been found in the literature. 

Sources of the various forms of interference may be studied by 

(1) Varying the angle (p 

(2) Utilizing both light beams or the test beam only 

(3) Varying the object plane of focus 

(4) Utilizing a monochromatic or a white (incandescent) light source 

(5) Varying the lateral extent of the light source 

(6) Varying the orientation and pitch of the Zehnder-Mach inter- 
ference fringes 

Figures 5 and 6 were obtained with monochromatic light at grazing 
incidence to an optically flat aluminized front-surface mirror 7 inches 
in length. Additional experimental conditions associated with the 
photographs were as f ollows : 

(1) Object plane was at the mirror midspan plane. 
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(2) Photographs in order a, b, . . . correspond to increasing 

cp, where photograph a corresponds to (p « 0. 

(3) Figure 5 involves only the test beam and figure 6 includes both 
the test beam and the reference beam. 

The wire appearing in the figures served to locate the surface at 
the object plane, and the irregularities (insulation) along the wire 
served to permit more critical focussing. In figure 5 the weak vertical 
bands were caused by glass striations in the camera lenses and are ir- 
relevant to the present discussion. 

In figure 6(a), for (p « 0, narrow reflection fringes similar to 
Lloyd's mirror fringes are visible. They eire produced by the interac- 
tion of the reference beam and the reflected beam (rather than the test 
beam and the reflected beam, which corresponds to Lloyd's Interference) 
because the fringes are inclined to the surface. As (p is increased 
(<p < 0), a second pattern of fringes (fig. 6(b)) gradually replaces 
the first pattern. The initial stage of formation of the second pat- 
tern is apparent in figure 6(a) . The second pattern is inclined to the 
surface at a considerably greater angle than the first pattern. The 
second pattern is probably caused by interaction of the first pattern 
with light diffracted by the plate edge nearest the camera because the 
fringes are curved (rather than straight, as would be expected if in- 
terference were the sole cause) . As <p is further increased until 
<p > 0, the second pattern disappears and a third pattern appears in- 
clined oppositely to the second pattern (figs. 6(d) and (e)). The 
third pattern is apparently caused by the interaction of light dif- 
fracted by the plate edge nearest the light source with the reference 
beam. The fringes are diagonal and curved and are absent when the 
reference beam is eliminated. When (p = 0, the spurious pattern, ex- 
cept for the diffraction pattern produced by the plate edge nearest 
the light source, should disappear if the light source is a geometrical 
point source. In practice, the light source is always of finite lateral 
extent. Thus, the diagonal fringes do not disappear completely, but, 
rather, a "picket-fence" effect caused by the combination of the two 
opposed patterns of diagonal fringes is observed (fig. 6(c)) when the 
light beam and the siurface are alined. 
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APPENDIX D 


END COERECTIONS 


The free -stream density p within a wind tunnel can he detennined 
from interference measurements if the boundary layers formed by gas flow 
along the wind-tunnel windows are taken into account. 

The basic equations for optical-path difference A\lf and fringe 
shift N are^ respectively. 


where s is a physical length measiired along a ray trace, and t is 
an optical-path length. 

For the present, consider a coordinate z that is parallel to the 
interferometer optical axis and has its origin at the wlnd-t\mnel wall 
which is initially traversed by the light. Let z increase positively 
in the direction of light propagation, as shown in flgxore 7. Let a 
second coordinate C = z be associated with an individual ray trace. 
Also, let = t-‘-2(z) represent the refractive-index distribution in 
the boxmdary layers adjacent to the wind-t\innel windows. For a ray 
trace which traverses the free stream and the window boundary layers. 


where 6„ is the window boundary- layer thickness, and the flow is as- 

simied to be symmetrical about the wind-tunnel midspan plane. For a 
ray trace through the reference hole r, which is connected to the ex- 
ternal atmosphere. 



(Dl) 


NX = 


(D2) 



= ^^r^ 


Therefore, 


N „X = i|r - \|f = p (L - 25 ) + 2 

^,r «, r «e ^ 


/: 


6 , 


z 
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Boundary layers formed on wind-tunnel windows are usually turbiilent. 
Thus^ assume that m- 2 is of the form 

^ \l/n 

M-g = l-^(O) + [m«> ■ h(0) j In the vicinity of 2=0 


Hg = fi(L) + [(I. - ^(l) j l\ 


l/n 


in the vicinity of z = L 


r^z r^z 

pgd^ = 

^^2,w + - t"2,w) 

o 

o 

L \ V J 


where p.(y = O) = d(0); eind^ because of symmetry, P-(o) = |a(L). Then, 

~ ^2,w^z ■ ^2,w^ n + 1 

where 

>^2,w =1^(0) = ^^(L) 

It follows that 


K.r'^= 


- 28 J - + as, 


- - l^r) (l- (>*2,W - >*r) ^ h 


However, 


|i = 1 + kp 


(D3) 


and the. optical-path difference in hole w with respect to hole r is, 
according to equations (D1) and (D2) , 


^w,r^ ~ 


Therefore , 


p - p = N JL 
w r w,r kL 


(D4) 


r 
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Assme that the densities at the windows and the other tunnel walls or 
model surface are equals that is^ that P2 w ” *^w* substituting 

expressions (D3) and (D4) in the previous equations for and 

solving for p yield 

•0 


p_ = 


p + 


W 


kL^ 


n + 1 


N, 


w,r 


kT.T.^ 


(D5) 


where 
= L - 


n + 1 z 
The method for determining p 


and p^ is as follows; 


(1) "White-light" interference fringes are oriented perpendlciolar 
to the wind-tunnel wall or model sin-face. Flow is then initiated in the 
wind tunnel. When equilibrium conditions are attained, the pressin-e, 
temperature, and density within holes w and r are, respectively, 
denoted by 


(a) Static hole w: Pw^ Tw? Pw 

(b) Reference hole r: p^, T^, 

where the wall or model temperature is assumed to be constant in the re- 
gion of interest. 

(2) The fringe pitch is adjusted until the zero-order fringe ap- 
pears centered in both holes simultaneously. This is always possible 
except when the flow is transient, as occurs in blow-down tunnels. If 
p^ = p^, the infinite-fringe adjustment results. 

(3) Alternate interferograms are obtained by using the desired 
monochromatic light and a white-light source. The white-light inter- 
ferograms are used for locating the zero-order fringes in the inter- 
ferograms obtained by using the monochromatic source when minor flow 
instabilities exist. Simiiltaneously, measurements of Pa^ Ta^ and T^- 
are made by using a barometer, a thermometer, and a thermocouple in- 
stalled at the surface in the region of interest, respectively. 

The values of the densities p and p can be computed in the 

m 

following manner: 




3378 


3378 


NACA TN 3507 


49 


r^ 

I 




(l) The atmospheric density is calculated from 

_ Pa_ 

Pa “ RT„ 
a 


where R is a gas constant^ and and Tg^ axe given by the barometer 

and thermometer, respectively. 

(2) Then, p^, is given by 


Pr 



Pa 


where T^ is obtained from the thermocouple. 

(3) The density p^ at the surface of the model is then calculated 

from 


Pw 


P-r + N, 


w,r 


_X^ 

kL 


where the number of fringes in the x' interval between the 

zero-order fringes in holes w and r, may be determined readily be- 
cause the fringe pitch d' may be measured and is independent of the 
magnitude of p in the two holes. 


( 4 ) The free-stream density p is given by 


P_ = 


Pp + 


N 


kL 


N, 


X5„ 


n + 1 w,r kT.T. 


where 



2 

n + 1 


5 

z 


and is the number of fringes in the x' -interval between the 

zero- order fringe in the reference hole and the zero- order fringe in 
the region corresponding to the free stream. 
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APPENDIX E 

REFRACTION AND CORNER CORRECTIONS (MODEL SPANS WIND TUNNEL) 

Evaluation equations which include corner^ as well as refraction, 
corrections can he derived hy extending the analysis contained in ref- t 
erence 1. 

Assvime that Cartesian coordinates x,y,z form a right-handed 
system having its origin at the intersection of the leading edge of the 
model with the wind-t\innel w a ll which is initially traversed hy the in- 
terferometer light heam. Coordinate z is taken essentially parallel 
to the light path, as shown in figure 8, and coordinates x and y 
are essentially perpendicular to the light path. A system of conjugate 
coordinates x',y',z' are determined in the camera- image space hy the 
relations 

x' = -mx 
y' = -my 
z ' = m^z 

where m is the lateral magnification produced hy the camera lens . A 
third set of Cartesian coordinates is associated with each in- 

dividual ray trace such that 



where the initial coordinates of any given ray trace entering the wind 
tunnel are (x,y,z) = (xQ,yo,0), or (4 ,ti,U = (0,0,0) . Consider houndary- 

layer refractive-index distrihutions : 

M- = l-i-Cy) adjacent to the model s\irface 

= ^ 2 ( 2 ) adjacent to the wind-tunn'el walls 

Let 5„ and 6„ denote the boundary- layer thicknesses associated with 
y ^ 

p and \i 2 , respectively. The quantities 5y and will always he 

regarded as positive. Although the geometry of the corner boundary 
layers may be similar to the geometries shown in reference 22, for 
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simplicity of the subsequent analysis the isopycnic lines will be assumed 

to form right angles in the corners. Thus, p and P2 extend into the 

comers and intersect along the diagonal siirfaces 

y/z = (0 < z 5 6^) 

y/(L - z) = by/b^ (L - z ^ z ^ L) 

Note that the fillet adjacent to the aunbient region is thereby neglected. 

The preceding assumption shoiild be reasonably valid because Kennard 
(ref. 4) has shown by example that a measured distribution p(y) is 
relatively insensitive to variations of the assumed corner geometry. 

The resiilting geometry and representive ray traces are shown in fig- 
ure 8. As in appendix symmetry of the boundary layers is assumed. 

As in reference 1^ p is assimied to be monotonic and representable 
for each ray trace in terms of j\ as 


= Po + apt^ H = £ ^v|y - yoT 

v=0 V =0 


in an interval y^ - where y = y^ when t] = p = Pq at t] == 0^ 
and ZV-yi <<1 because p « 1. The subscript notation is used for t] 

in order to differentiate between p as a function of y and p as a 
function of t]. The coefficient b^ is defined by 


b 


V 




where the choice of the + or - sign depends upon whether p is an 
increasing or decreasing fimction of y^ respectively. 


Light Path 

The ray traces are straight lines ^ except in the field p^ where ^ 
as in reference 1^ it is assumed that each ray trace can be represented 
by the power series 


«• 

E 

a=0 



a 
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where ^ hy definition^ 





If the light source is an axial-point source^ the limiting ray-trace 
displacement t] = H in the field p is given by 



* cr=2^4^ 6^ • . • 


where 


r.* ^ -k 

^2 = 2 


♦ 1 V 

- 12 ^ 1 ^ 2 ^ • ' • 

and the reduced model span Lp is given by 


Is - L - (El + Eg) 


where, by virtue of the boundary conditions, 

6 

u z 

5i = ST ^0 

<J 

^2 “ 6y 

The ultimate displacement at C = L is, approximately. 


Distortion 

In closed form the distortion is given by 


D = H + (?2 - 
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where 

K = 1 - ^ 

L 

and zp is the real locus of points which are imaged at the selected 

image plane. By virtue of the ray- trace series and the corresponding 
expansion 



0=2 ^ 4 ^ 6 ^... 

for its derivative^ the distortion is given in series form by 

- a(KL - 

0 = 2 ^ 4 ^ 0 ^ . . • 



Fringe Shift 

In the present instance the measured fringe shift serves to as- 
sociate the boundary-layer density distribution p(y) with the free- 
stream density p . By applying eq\iations (D1) and (D2) , the optical 

length of a ray trace through the plate and window boundary layers 

is given by 




r v^i 

■' L 

^lfl = 

+ 

p ds + 

P2<is 

Vj 

0 

V. 

Cl J 



whereas for a ray trace through the free stream and window boundary 
layers , 

- 25^) + 2 + p^KL(sec 3^ - l) 


where 



tan 3 = d'o/d-C 


IT 
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P = Pj at ^ and the last tem in the expression for ^ results 

from refraction (ref. 20). Ray traces through the field \i are given 
in terms of the auxiliaiy coordinates by the differential 

equation (ref. 20) 




It follows that 


Thus^ NX =5 or 


sec p = 


_ll 


^0 


'^1 fV^l 

NX= I + 2 

'o 





- H (L - 2& ) - 2 


J 


f5 


0 


- iT ^ A4g 


0 


(El) 


which results by substituting the right sides of the four preceding ex- 
pressions into the expressions for \lr.j and and then expanding 

2 

1-1 /n^ in the form 
n' 0 


li. 


^0 


= p. + 2Ap + . . . 

0 n 


where terms in Ap^ . . . are infinitesimals of higher order . 

Algebraic expressions for those integrals which Involve pg depend 
upon the functional form of p„. However^ the integral involving Ap 

Tj 

can be evaluated in series form^ as in reference 20. 

Solution for an effective-average wlndow-bovindary-layer density 
distribution . - Assume that form 
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^I2(A)](B - A) 


which is strictly not an average density assumption although it, hence- 
forth, will he referred to as such. Because of symmetry, 
p(0) = p(l) = equation (El), therefore, becomes 

^a 

^ + l^2,w)^2 - - 25z) - (|i2,w + ^ KL Ap^ 

The fourth term on the right side of the preceding equation can be ex- 
panded as follows: 

2p^ W ^2,w^^2 ^ 2 (j^O + ^^2,w ~p^ ^^^2 

However, 

^2,w *^a 

77" * 1 


Thus, letting 

^ 2 ^^1 ■*■ ^2^ = ^ - 2 (^1 + Cg ) 

there results, following rearraingement of terms. 


NX - (Pq - - ^^2,w^ [\ ■ 2 ^^1 ^2^] 


2 



Ap^d^ + I (3^2 " 2KL)Apg 


By utilizing the series expansions for \i and and by translating 
the coordinates of the integral in order to give a lower limit of zero^ 
the preceding equation becomes in series form 
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NX 


= (‘‘o - ^ - I + ^2)] + 


3C E ^1^-. 


v=l ff=2^4,6, 


i(352-2KL)'^ C ^=a,.4 

v=l a=2,4,6^... 

where the coefficients ^a,v given hy (ref. 20) 

^2^1 == i *^2,2 = ® 


^4,1 =12 ^1^2 


1 1.2 

^4,2 - 4 ^1 


(E2) 


Solution for power-lav window-boimdary- layer density distribution , - 
Assume that 



. 'U 

(0<i< 5 ^) 

(L - 62 ^ C = L) 


Because of symmetry^ l^(0) = I-i(L) = t>efore. Following integra- 

tion of those integrals which involve 10.2 (in the same manner as in 

appendix D) and the collection^ whenever possible, of like terms in p, 
equation (El) becomes 


NX - (1^5-1^.) i Ci+?2-2^- HTT 


n+1 n+1 

•©■ - ! 




^1 


Ap^d^ + 


1^2, W ^ ^‘■a . II 

Pq ^2 n+1 Pq 


n+1 n 


Apij 


2 
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GO 

t 


However^ loecause 


.1 


>*0 “o 


and the term — 

adjacent terms ^ It follows that 

NX = (^Xo - + 


n +1 

n 


82 is of smaller order than the 


- "a,wM^TT \ * ^2) + 5^ \ 


n +1 n+l ~i 


Lr-Ki 


+ (C2 - KL)AHg 


which becomes in series form 
NX = (^ 1 q - ii^)Lp + 




n + 1 * ^ 2 ^ ■'■ n + 1 ^’z 




•n +1 n+l-i 

^^l\^ /^ 2 \“ 


E E 


v=l a=2^4_,6, . . . 


v=l o=2j4r^6j . . . 


he II 

V 


(E3) 


Evaluation Equations 

The final evaliiation equations can he derived hy applying the pre- 
ceding results in the same manner as in reference 1. The procedure is 
as follows; After 1 + kp is substituted for p (with the appropriate 
subscript retained), equation (E2), or (E3) , is solved for p* = p^/p , 

O' •D'^ 

the choice of equation depending upon the choice of the window-boundary- 
layer assumption. The quantity Pq represents the density associated 

with the profile p(y) at the ordinate value jq . The expansion 
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determines the final form of the evaluation eq.ua tions. Thus, the final 
equations are denoted according to the mmber of terms of the preceding 
series which are utilized. Only the one- or two-term approximations 
appear to he of practical use when corner corrections are to he made. 

By considering only two terms of the series, all coefficients 
h^ (v >l) are assumed to he zero; hence all coefficients c* and 

C(j^ V which involve h^ (v > l) are also zero. If only the first term 
ho = do of ttLe series expansion for is considered, then all the 

Z sums are zero. Because p* = at y = y^, whereas the corre- 

sponding measured value of N is associated with y = y^, the indi- 
cated distortion correction must he performed in order to determine the 
value of yo« If the preceding operations are performed, the evaluation 

equations (with the exception of h^) presented in the section Density 

Field p(y) (Model Spans Wind Tunnel) are obtained. The derivation of 
expressions for hj_ is considerably more complicated. First, equation 

(E2), or (E3) , is differentiated with respect to y. An expression for 
dW/dy at y = jq is, thereby, obtained. However, by measurements on 

an interferogram a distorted value dN/dyjj at y = y^ is obtained. 

It is shown in reference 1 that the two quantities are related by 


Therefore, the series expansion for D must also be differentiated 
with respect to y. When the expressions for dW/dy and do/dy are 
substituted in the preceding equation, three additional assumptions, 
namely. 



(1) b^L « 1 

(2) 62/^ of order of unity 

(3) biL62/5y « 1 


permit reduction of the eqiiation to a form which is linear with respect 
to bj. Because (dri/d^)^ = bj_L is the maximum slope of any given ray 
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trace with respect to the optical axis, assumption 1 is certainly valid 
in practice. Moreover, if the window-boundary-layer thickness does not 
exceed the model-bo\mdary- layer thickness by an order of magnitude, then 
assumptions 2 and 3 are satisfied. The intermediate equations are 


lb 


_ 

1 ^ dyp ^ 



+ kp 




for the effective-average density assmption, or 



for the power-law density assumption. In practice all terms in each 
expression for bj are usually small compared with the first term (say 

less than 10 percent of the first term) and, also, will tend to cancel 
when algebraic signs are considered. The final expressions for bq 
are then, respectively. 


b 


1 


dN ^ 


or 


bi = 


dN 


% 


In any given experiment the maximum values of the neglected terms should 
be estimated in order to confirm the validity of the final solutions. 
Unless the resultant magnitude of the neglected terms is considerably 
greater than 5 or 10 percent of the first term, they can still be neg- 
lected because the error in determining du/dyjj alone can cause an er- 
ror of 5 or 10 percent in determining b^^. 


j 
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APPENDIX F 


REFRACTION AND CORNER CORRECTIONS (MODEL NOT COIWAINED 

WITHIN WIND tunnel) 

If the model is not contained within a wind timnel^ then the func- 
tion the coordinate systems^ and the quantities 6y and 5^ can 

be defined as in appendix E. The additional fimction \i^ = M* 3(y^z) is 

introduced to represent the refractive- index distribution at the corners 
and defines the extent of the corner boundaiy layers in the 

z -direction. The appropriate geometry and representative ray traces 
are shown in figure 9. The assimied geometry of the comer boundary 
layers is likely to be more critical than in appendix E because the 
radii of the isopycnic lines around the ends of the model are relatively 
larger. However^ Kennard^s results (cf. appendix E) apply specifically 
to the present case. It will be assumed that the outline of the bound- 
ary layer in the vicinity of the spanwise ends of the model is repre- 
sentable by the hyperbolas 

(y - 5 y )(62 + z) = -a^/2 

in the vicinity of z = 0, and 

(y - - L - S^) = a^/2 

in the vicinity of z = L, where a is the perpendicular distance from 
the isopyncnlc line bounding the ambient fluid to the intersection of 
the asymptotes of the hyperbola. 


Light Path 


In the field p, 

® - E 

cj— 2 • . • 

where the reduced model span is given by 

L, . L - - S^) 
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The ultimate displacement in the space-variant field is, approximately. 



Distortion 

In closed form the distortion is given by 

D . E - (KL + 8, - 

/L^ 


or^ in series form^ 


® " * E [% - »(KL + 8^ - 

0=2,4, 6, . . . 


Fringe Shift 


The determination of the boundary- layer distribution p(y) with 
respect to the ambient density is of interest. Thus, 


'iz-5z 




= 


_ + 


''L+b^-Zz+^g 


U ds + 


■^z“^l 


^z-&z 




Pjds 


{(L - 2 Z 2 ; + 282 + + C 2 ) (Kb + ^2 “ Z2)(sec - 1^ 

The upper limit of /^^ds is only approximate • However^ because 

(1) The exact profile of the boundary-layer outline is indefinite 

(2) The exact corner density distribution is unknown 

the assumed upper limit appears to be sufficiently accurate for practi- 
cal purposes. Performing the operations indicated in appendix E yields 
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NX= = ^^0^ + 2 


L-L, 


An ds + 


^ L-L. 


L-Le-Ci 


n^d^ + — 


Le+Cs 




>*.(% + 5l + £2) - ^^[(K - 1)L + I-e + £2]^1': 


H 


(Fl) 


where 


== L - 


7+6 
z z 


Assuming an effective-average end density distribution corresponds 
to assxxming that 


L-Lg 

= i (n^ + n^)I 




^e^2 


= 2 + O^z 


Then equation (Fl) reduces to 


^ L 


NX = (Hq - n_)Lj, + 2 


_ A4 « -f(K - 1)L + - i 

L-L^ 


where the effective- test- section span is 


Lj, = Is * I (£1 y 


In series form the result is 


NX = (^^0 - 



V =1 


i : 

cj— • • • 


\ -^+1 
a + 1 ^ 
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Evaluation Eq.uations 


The final evaluation equations are obtained in the same manner as 
in appendix E. If, as in appendix E, it is assumed that biL « 1, 
then 



L(yo - Sy)'^ 


(y- 


H 




which reduces to 


b 


1 


dU 


when the second term is small compared with the first term. The second 
term is usually small compared with the first term because, by assiamp- 
a2 

■tion, — —j — — < I This assumption is required by the 

2(yo - 

fact that the isopycnic lines are assumed to be hyperbolic at the' 
corners. It can only be violated when yQ, yjj -*■ 6y However, in the 

vicinity of y = 5y, p(y) -*■ and the entire corner correction is 
then known to be negligible from physical considerations. 
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APPENDIX G 


EFFECT OF WONUWIFORM HEATING OF WIND-TUNNEL WINDOWS 

Consider two ray traces which traverse a test section hounded by 
windows. Assume that the traces are parallel to the optical, or z-, 
axis and that refraction is negligible. Let t(z) represent the tem- 
perature difference at any given value of z along the two paths 
within each window, where by symmetry t(z) possesses the same dis- 
tribution throughout both windows. Then, neglecting possible ef- 
fects of thermal stresses within the glass , 

At = 2(apg + p) T(z)dz 

within the two windows, where t is the window thickness. The co- 
efficient of thermal expansion a and the temperature coefficient of 
refractive index p possess, for ordinary crown glass (Pg « 1.52), the 

respective values 


a * 8X10"®/°C 
P * 1.5X10“®/°C 

at 20° C for X = 5893 A (sodiimi D-lines) . In terms of the error Ap* 
in p*, the preceding expression becomes 

2(c4» + 13) P 

Jo 

Usually “tCz) « |T2,w - so that t(z) can be regarded as a linear 

function of z, corresponding to heat conduction through a plane wall 
of Infinite lateral extent. Thus, 


Ap* = 


2 (op. 


kP. 


ilp(O) 


where x(o) represents the temperature difference at the internal surfaces 
of the windows. The preceding values of the required constants were used 
to compute the value of Ap* reported in the section Supersonic Air Flow 
Along a Flat Plate . 
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APPENDIX H 


CYLINDRICAL MODEL 

The cylindrical model consists of a 3-inch-dlameter , 8-inch-long 
Duraliamin tube with 1- inch- thick walls and is bounded on both ends by 
1-inch-thick removable end caps. The end caps support a 3/4-inch- 
diameter alundum tube woiond with nichrome resistance wire within the 
outer cylinder. The rests tauce-wire leads were passed out through 
the centers of the end caps (one lead through each cap) and insulated 
from the caps by Fiberglas sleeving. The spacing of the resistance 
wire was reduced near the ends of the tube in an attempt to produce 
uniform heating of the exterior siarface of the outer cylinder. Iron- 
constantan thermocouples were located at 11 points on the exterior sur- 
face of the Duralumin cylinder. Temperatures indicated by the thermo- 
couples were recorded automatically. 

For experimental piorposes the heater was suspended from an ad- 
justable horizontal supporting rod by piano wire. Thus^ the heater 
was, effectively, freely suspended in space. Air conditioning in the 
large room in which the experiment was performed was shut off, and the 
cylinder was shielded from stray air currents. The heater current was 
controllable by means of an autotransformer, thus permitting attain- 
ment of a virtual steady- state condition while obtaining an Interfero- 
gram and simultaneously recording the surface temperatures . 

A small ambient temperature rise (approximately 1° F) was recorded 
in the shielded region during each run by means of a thermometer. The 
effect of the ambient temperature rise was included in the calculation 
of the temperature profiles. The temperature variation over the ex- 
terior surface of the cylinder was found to be less than 2° F through- 
out the range 24° F < 0 < 450° F. 
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(a) Orientation of fringes parallel to grad p . 



(b) Orientation of fringes perpendicular to grad p . 
Figure 1. - Interference-fringe orientation. 
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(a) Standard deviation of a single measurement of pitch of constant-pitch Interference fringes 
as a function of mean Interference-fringe pitch. 



0 .4 .8 1.2 1.6 2.0 2.4 2.8 3.2 3.6 4.0 

Mean Interference-fringe pitch, , mm 


(b) Standard deviation of a single measurement of unit fringe shift from measurements ">n constant- 
pitch Interference fringes as function of mean Interference-fringe pitch, 

Figure 2. - Error In measuring fringe shift. 




Fringe-Shift derivative, dN/dyp, In.-l Fringe shift 
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(b) Comparison of hypothetical (ref. l) and optically measured profiles dNjy'dyp. 
Figure 3. - Comparison of hypothetical and measured profiles dNj>/dyj^. 
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Test beam and 
reference beam 



(a) cp < 0. 


r 


I 


Test beam and 
reference beam 
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Diffracting edge^j^ Model ^ 


Reference 
beam only 


(b) 9=0. 


z 



Test beam and 
reference beam 


Reference 
beam only 


Figure 4. - Ray-trace orientations. 
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(c) <P >>0. (b) <P=*0. 

Figure 5. - Surface effects, test beam only. Fraction of test- section span K, 1/2; magnification, X28. 


(Lm 


L 


MCA TW 3507 








(c) <P=0- 



(c) <p=0- 


C-38496 


Figure 6. - Spxirious fringes, both beams. Fraction of tost- section span K, l/2; magnification, X19. 
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T = T = T„ 

P = Pr P = P;, 

(b) Longitudinal section. 


Figure 7. - Wind-tunnel geometry for determining end corrections. 
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Figure 9. - Kay traces in proximity to model not contained in wind tunnel. 
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Figure 10. - NACA 6- inch Zehnder-Mach interferometer vith schlieren provision and afocal 
camera-lens system. 
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Figure 11. - NACA 6- Inch Zehnder-Mach interferometer (side plate removed). 
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(a) Density profiles; x * 2.5 inches. 

Figure 12. - Measurements of laminar boundary layer associated with supersonic flow along an insulated flat 
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(b) Velocity profiles J x = 2,5 inches. 

Figure 12. - Continued. Measurements of laminar toundary layer associated with super- 
sonic flow along an insulated flat plate. 
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(c) Section skin-friction coefficient as function of section Reynolds number. 

Figure 12. - Concluded. Measurements of laminar boundary layer associated with 
supersonic flow along an insulated flat plate. 
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Figure 13. - Comparison of successive approximations of density profile in flat-plate boundary layer 
(supersonic flow). ~ 0.989x10®. 
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Figure 14. - Boundary- layer cross sections. 
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Figure 15. - Details of flat-plate model used in subsonic -flow 
experiment (ref. 24) . 
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(a) Unfiltered mercury light. 



(b) Filtered mercury light (X. = 5461 S.) . 

Figure 16. - Typical interferograms of boundary layer induced by 
subsonic flow along flat plate. 
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Figure 17. - 
profiles . 


Subsonic air flow along a flat plate. 
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Figure 18. - Orientations of light heam with respect to cylinder. 
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(a) Orientation 1: Light propagation parallel to cylinder axis. 

Figure 19. - Typical interferograms of heated horizontal cylinder. 
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(b) Orientation 2: Light propagation perpendicular to cylinder (c) Orientation 3: Light propagation perpendicular to cylinder 

axis at mldlen^h plane. axis at end plane. 

Figure 19. - Concluded. Typical Interferograms of heated horizontal cylinder. 
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(a) Two-term approximation including corner corrections; K = 1/2; 

= 83.3° F; v based on . 

Figure 20. - Dimensionless temperature profiles beneath horizontal 
circular cylinder. 
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(b) Wachtell-DePrate method (refs. 2 and 3); K = l/3; = 80.2^ F; 

V based on T^. 

Figure 20. - Continued. Dimensionless temperature profiles beneath 
horizontal circular cylinder. 
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(cj_ Two-term approximation Including corner corrections; K = 1/3; 

=» 80.2° F; v based- on T^. 

Figure 20. - Continued. Dimensionless temperature profiles beneath 
horizontal circular cylinder. 
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(d) Van Voorhls-Wachtell method (refs. 29 and 2); = 79.9° F; 

V based on . 


Figure 20. ^ Concluded. Dimensionless temperature profiles beneath 
horizontal circular cylinder. 


Temperature ratio, 


94 


MCA TW 3507 



Figure 21. - Dimensionless temperature profiles beneath horizontal 
circular cylinder. K = 1/2; T„ = 83.3° F; v based on T^. 
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Figure 22. - Interferogram of hot-plate model described 
in reference 1; K=l. 
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